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FOREWORD 


This  is  the  fourth  Interim  Technical  Report  on  Contract  F33615-76-C-1198, 
covering  the  contract  period  of  1  May  1978  to  1  April  1979.  The  first  interim 
report,  AFAL-TR- 77-17,  covers  the  period  4  May  1976  to  4  November  1976;  the 
second  interim  report,  AFAL-TR- 77- 198,  covers  the  period  4  November  1976  to 
4  May  1977;  the  third  interim  report,  AFAL-TR- 78- 31,  covers  the  period  of 
4  May  1977  to  4  November  1977. 

The  reported  work  was  performed  in  the  Bubble  Memory  Technology  Department, 
Electronics  Research  Center,  Rockwell  International  Corporation,  Anaheim,  CA. 

Dr.  G.  R.  Pulliam  is  Department  Director  and  Program  Manager  of  this  contract. 

D.  M.  Heinz  is  Principal  Investigator.  V.  S.  Speriosu  is  the  major  contributor 
to  this  report. 

On  1  May  1978,  a  study  was  initiated  on  the  physics  of  ion- implantation  in 
magnetic  garnets  with  the  aim  of  determining  requirements  for  simple  propagation 
structures  in  contiguous  disk  types  of  bubble  memory  devices.  After  work  on  this 
contract  task  ended,  the  investigation  was  continued  on  a  Rockwell  International 
internally- funded  program.  During  this  period,  the  use  of  ferromagnetic  resonance 
measurements  were  developed  to  complement  the  X-ray  diffraction  techniques  evolved 
earlier.  The  combined  results  of  the  contract- supported  and  Rockwell  supported 
studies  are  presented  in  this  document  to  provide  a  complete  account  of  this 
investigation. 

The  contractor’s  report  number  is  C76-845. 47/501. 
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SECTION  I 


INTRODUCTION 


The  objective  of  this  research  effort  is  to  demonstrate  the  feasibility 
of  a  new  technology  for  all  solid  state  large  data  base  memories  for  airborne/ 
spacebome  applications.  The  program  scope  is  limited  to  demonstrating  the 
feasibility  of  the  device  functions  and  bit  density  capabilities  required  to 
meet  the  system  goals. 

It  is  the  objective  of  this  program  to  advance  one  technology  into  de¬ 
velopment  .  Based  on  an  evaluation  of  the  present  and  potential  capabilities 
of  candidate  solid  state  technologies  to  meet  the  system  goals,  magnetic 
bubble  domain  technology  has  been  selected  as  the  one  to  be  pursued  on  this 
program. 

As  the  bubble  device  period  of  a  permalloy  propagation  structure  is 
scaled  down  to  attain  higher  bit  densities,  two  serious  problems  arise:  one 
is  related  to  photolithographic  resolution  and  the  other  is  related  to  in¬ 
plane  magnetic  drive  field.  Currently,  the  most  advanced  permalloy  propaga¬ 
tion  structure  is  the  asymmetric  chevron  whose  minimum  feature  is  about  one 
eighth  of  the  circuit  period.  Since  the  limit  of  optical  photolithographic 
resolution  is  >lvim,  a  4ym  period  device  must  be  fabricated  by  some  other 
lithographic  technique  such  as  electron  beam  microfabrication.  The  minimum 
drive  field  of  a  permalloy  device  is  proportional  to  the  magnetization  of 
the  bubble  material  (which  is  approximately  inversely  proportional  to  the 
bubble  diameter).  As  a  result,  the  drive  field  of  a  permalloy  device 
increases  rapidly  with  decreasing  bubble  diameter.  Existing  8pm  poriod 
permalloy  devices,  though  not  optimized,  require  more  than  60  Oe  operating 
drive  field.  Thus,  even  when  the  electron  beam  microfabrication  technique 
becomes  available,  the  expected  increase  in  drive  field  casts  some  doubt  on 
the  viability  of  4pm  period  permalloy  devices. 

The  ion- implanted  contiguous  disk  device  appoars  to  solve  both  of  the 
above-mentioned  problems.  As  far  as  the  propagation  structure  is  concerned, 
the  minimum  feature  of  an  ion- implanted  device  can  be  increased  by  a  factor 
of  three  (or  more)  over  that  of  the  asyirmetric  chevron,  which  makes  possible 
the  fabrication  of  a  4pm  period  device  by  conventional  photolithography. 

In  contrast  to  the  permalloy  device,  the  ion- implanted  device  has  no  circuit 
gaps  along  its  propagation  track.  Consequently,  the  energy  of  a  bubble  is 
nearly  independent  of  its  position  along  the  propagation  track,  resulting  in 
a  very  weak  dependence  of  the  minimum  drive  field  on  the  garnet  magnetization. 
Thus,  both  the  lithographic  resolution  and  the  high  drive  field  problems  of 
4 vim  period  permalloy  devices  may  be  solved  by  ion- implanted  devices. 
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Ion- implantation  is  an  established  technique  for  changing  some  of  the 
properties  of  as-grown  crystals,  having  been  successfully  applied  to  both 
semiconductors  and  garnets  (Ref  1).  In  recent  years,  ion- implantation  has 
assumed  a  key  role  in  the  treatment  of  magnetic  bubble  materials  for  hard 
bubble  suppression  (Ref  2) ,  as  well  as  in  the  development  of  contiguous  disk 
bubble  memory  devices  (Ref  3) .  Although  widely  used,  the  effects  of  ion- 
implantation  are  not  very  well  understood.  No  a  priori  guidelines  for 
selecting  optimum  implantation  conditions  are  available.  For  contiguous 
disk  devices,  the  final  test  is  device  performance,  characterized  by  minimum 
drive  field  and  maximum  bias  margin.  Since  the  trial -and -error  approach  for 
optimizing  implantation  by  preparing  processed  devices  would  be  prohibitively 
time  consuming  and  expensive,  the  more  fundamental  approach  of  investigating 
the  physics  of  ion- implanted  garnet  was  adopted. 

A  study  was  initiated  on  this  program  to  determine  the  dependence  of  the 
garnet’s  crystalline  and  magnetic  properties  on  ion- implantation  conditions: 
Species,  energy,  dose  and  temperature.  One  of  the  major  macroscopic  changes 
in  the  crystal  structure  is  the  expansion  of  the  unit  cell,  which  places  the 
implanted  garnet  in  compressive  stress.  In  bubble  materials  with  a  negative 
magnetostriction  constant,  the  stress  induces  a  negative  magnetic  anisotropy 
field,  Hi,,  whose  value  and  depth  dependence  greatly  influence  device  perform¬ 
ance.  The  microscopic  displacement  of  atoms  from  their  perfect -crystal 
positions,  usually  termed  implantation  damage,  affects  other  magnetic  pro¬ 
perties  such  as  magnetization,  4rrMs,  coercivity,  Hc,  exchange  constant.  A,  and 
gyromagnetic  ratio,  y .  Knowledge  and  control  of  these  properties  and  of 
their  stability  over  a  wide  temperature  range  are  essential  to  device  manu¬ 
facture  . 

Several  methods  are  presently  in  use  for  the  study  of  ion- implantation: 
The  rate  of  chemical  etching  is  known  to  depend  on  the  extent  of  damage 
(Ref  4) ,  the  presence  of  additional  absorption  peaks  in  spin-wave  resonance 
is  interpreted  as  an  indication  of  a  layer  with  negative  anisotropy  (Refs  5, 
6) ,  changes  in  magnetic  bubble  collapse  fields  are  used  as  a  rough  measure 
of  the  implantation- induced  flux-keeper  effect  (Ref  7),  and  both  the 
Rutherford- scattered  ion  flux  (Ref  8)  and  low  energy  electron  diffraction 
(Ref  9)  are  related  to  lattice  structural  changes.  Although  each  method  has 
proved  useful  in  understanding  implanted  crystals,  their  common  feature  is 
the  difficulty  in  obtaining  detailed  and  quantitative  interpretations  of 
experimental  results . 

Double  crystal  X-ray  diffraction  has  also  been  applied  to  the  study  of 
ion- implanted  materials  (Refs  10,11,12).  Diffraction  profiles  (rocking 
curves)  contain  information  about  the  small  changes  in  lattice  parameters 
(strain) ,  microscopic  displacement  (damage) ,  and  other  inhomogeneities 
caused  by  implantation.  A  method  of  extracting  this  information  from  experi¬ 
mental  rocking  curves  has  been  developed  on  this  program  and  has  been  used 
successfully  on  implanted  garnet ,  as  described  in  later  Sections  (Ref  13) . 
The  method  is  general  and  applicable  not  only  to  implanted  crystals,  but  also 


to  diffused  or  as -grown  crystals  whose  crystallographic  properties  vary  with 
depth.  The  relative  simplicity  of  the  experimental  technique  and  mathematical 
treatment  should  make  this  approach  accessible  to  a  wide  variety  of  workers. 

In  addition  to  X-ray  diffraction,  ferromagnetic  resonance  (FMR)  has  been 
used  to  obtain  direct  values  for  the  anisotropy  field,  H^,  the  gyromagnetic 
ratio,  y,  the  damping  parameter,  a,  the  exchange  constant,  A,  and  rough 
values  for  the  magnetization,  4ttMs.  A  method  yielding  the  death  dependence 
of  in  an  ion- implanted  film  has  been  published  recently ,  which 
showed  that  the  profile  of  Hv  is  proportional  to  the  strain  distribution,  as 
*  predicted  by  theory.  This  FMR  method  has  been  adapted  as  a  complementary 

tool  to  our  X-ray  technique  in  analyzing  ion- implanted  garnets. 

The  remainder  of  this  report  consists  of  the  following  parts:  Section  II 
describes  the  kinematical  model  of  X-ray  diffraction  and  its  application  to 
implanted  garnet;  Section  III  presents  the  results  and  discussion  of  the  study 
of  Ne+,  He+,  and  If1-  implantations;  and  Section  IV  summarizes  the  results  of 
this  investigation. 


SECTION  II 


KINEMATICAL  X-RAY  DIFFRACTION  IN  THIN  CRYSTALS 


X-ray  diffraction  is  a  standard  laboratory  tool  for  the  determination 
of  crystal  structure.  In  the  usual  double-crystal  experimental  set-up, 
the  X-ray  output  of  a  metal  target  is  collimated,  monochromatized  and 
partially  polarized  by  a  first  crystal  of  good  quality.  The  second  (sample) 
i  crystal,  whose  lattice  parameters  are  to  be  measured,  is  slowly  rotated 

!  with  respect  to  the  incident  X-ray  beam,  while  the  diffracted  beam  intensity 

is  monitored.  The  intensity  profile  as  a  function  of  angle  is  called  a 
rocking  curve.  Figure  1(a)  schematically  shows  the  rocking  curve  of  a  thick 
crystal  with  uniform  lattice  constant.  The  intensity  is  zero  except  for  the 
narrow  angular  range  which  satisfied  the  Bragg  condition.  The  position  of 
the  peak  is  determined  by  the  lattice  parameter  of  the  as -grown  crystal, 

(in  this  case,  an  epitaxial  garnet  film  grown  by  liquid  phase  epitaxy,  LPE). 
The  width  and  detailed  shape  of  the  peak  indicate  the  degree  of  uniformity 
and  perfection  of  the  crystal. 

After  ion- implantation  (by  appropriate  implantation  species,  energies, 
and  dosages) ,  an  increase  in  diffracted  beam  intensity  is  observed  on  the 
low  angle  side  of  the  LPE  peak.  As  shown  in  Figure  1(b),  this  intensity 
has  an  oscillatory  character  and  is  typically  10  to  100  times  less  than  the 
LPE  peak.  The  period  of  oscillation  varies,  but  is  always  greater  than  the 
width  of  the  LPE  peak.  Since  the  thickness  of  an  implanted  layer  is  small 
compared  to  both  extinction  and  absorption  lengths (Ref  1S)>  a  kinematical 
model  of  diffraction  is  valid™e*  .  In  this  section  we  show  how  a 
kinematical  treatment  of  the  diffraction  problem  can  account  for  the  de¬ 
tailed  shape  and  intensity  of  the  implanted  layer  rocking  curve.  By 
adopting  this  approach,  quantitative  values  for  the  strain  and  damage  dis¬ 
tributions  may  be  obtained,  and  consequently,  the  thickness  of  ion- implanted 
layer. 

Kinematical  theory  deals  with  weakly  reflecting  planes  of  spacing,  d, 
(which  is  not  necessarily  constant  throughout  the  crystal) .  The  calculation 
of  the  total  reflected  amplitude  consists  of  the  addition  of  the  plane  wave 
amplitudes,  reflected  (only  once)  by  each  plane,  with  proper  account  taken 
of  the  phases.  Simplifications  result  if  the  planes  are  parallel  and 
infinite  in  extent.  Disorder  in  the  crystal  lattice  is  introduced  in  a 
manner  similar  to  the  Debye-Waller  treatment  of  thermal  fluctuations,  i.e., 

I  by  decreasing  the  reflectivity  of  each  plane.  This  is  equivalent  to  a  de¬ 

crease  in  the  structure  factor  of  a  given  reflection. 

The  strain  and  damage  distributions  in  ion- implanted  crystals  have 
,  been  found  to  vary  only  with  depth,  being  constant  and  isotropic  at  a  given 

'  distance  from  the  surface.  At  any  depth,  the  strain  perpendicular  to  the 

surface  is  defined  as  the  fractional  change  in  the  spacing  between  planes 
parallel  to  the  surface,  (d-d0)/d0,  where  d  is  the  strained  spacing  and 
is  the  undisturbed  spacing.  Except  for  very  high  implantation  dosages, 
when  amorphous  layers  may  be  created,  the  constraints  imposed  by  the  un¬ 
affected,  bulk  (film  and  substrate)  crystal  cause  the  strain  parallel  to 
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between  the  largest  implanted  layer  peak 
and  the  LPE  peak  is  about  0.3  degrees.) 


the  surface  to  be  much  less  than  that  perpendicular  to  the  surface. 


The  effects  of  ion- implantation  are  represented  in  our  model  by 
an  arbitrary  number,  N,  of  laminae  parallel  to  the  surface. 

Hach  lamina,  j,  has  its  own  number  of  planes,  n j ,  inter-planar  spacing,  d j , 
and  damage  function,  Wj .  The  damage  function  varies  from  zero  (amorphousJ 
layer)  to  one  (perfect  crystal) .  The  condition  for  the  validity  of  a  kine- 
matical  approach  is  that  the  thickness  of  any  lamina  i,  n^dj,  be  much  less 
than  L0,  the  extinction  length.  If  the  total  thickness  is  much  less  than 
the  absorption  length,  £0,  absorption  may  be  safely  neglected.  For  garnets, 
typical  values  of  L0  and  £0  in  strong  reflection  are  2ym  and  lOym,  respective 
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For  unit  incident  amplitude  and  a-polarization,  the  complex  amplitude 
reflected  by  an  undamaged  plane (Ref  17)  is  given  by 
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where  dQ  =  undisturbed  inter-planar  spacing;  X  =  X-ray  wavelength;  V  =  volume 
of  unit  cell;  0g  =  Bragg  angle  of  undisturbed  crstal;  e^/mc^  =  classical 


electron  radius;  F 


hkd 


structure  factor  for  reflection  (hk£) . 


The  total  amplitude  reflected  by  the  N  laminae,  obtained  by  summing  the 
amplitudes  coming  from  individual  planes,  account  being  taken  of  phase  differ 
ences  and  damage,  is  given  by 
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In  symmetric  reflections  and  for  ©b<tt/4 ,  the  amplitude  at  the  surface 
of  a  thick  perfect  crystal  (derives  from  Ref  15,  Eq.  35.1)  is 
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and 

c  = 

1  for  a-polarization  or  cos  20  for  r-polarization. 

Experimentally,  it  is  found  that  the  Bragg  peak  reflected  by  the  bulk 
LPE  crystal  is  broader  and  less  intense  than  the  predictions  of  the  perfect 
crystal  dynamical  theory.  This  is  probably  due  to  composition  and  lattice 
parameter  variations.  In  the  angular  range  of  interest,  corresponding  to 
the  tail  of  the  Bragg  peak,  this  discrepancy  can  be  eliminated  by  fitting 
the  measured  bulk  intensity  with  a  function,  such  as  an  exponential ,  and 
using  its  square  root  as  the  bulk  amplitude.  A  good  fit  can  usually  be 
obtained  by  simply  multiplying  the  theoretical  perfect  crystal  amplitude  by 
a  fitting  factor,  f,  whose  value  (between  1  and  3)  is  chosen  for  the  best 
agreement.  This  procedure  is  legitimate  since  we  are  not  interested  in  a 
theoretical  model  for  the  bulk. 

Neglecting  multiple  reflections  between  the  ion- implanted  layer  and 
the  bulk,  for  unit  incident  amplitude  and  cr-polarization,  the  total  reflected 
amplitude  becomes 

R  (e)  =  An  (e)  +  f  X(e)  e‘i6l  (4) 


and  the  intensity  is 

I  (e)  =  |R  (e)|2 


*  |An  (e)  | 2  +  [f  X(e)]2  +  2f  XU)  cos  \  .  (5) 
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Since  the  incident  beam  is  not  perfectly  o-polarized  by  its  reflection 
from  the  first  crystal,  the  top  layer  amplitude,  Am,  should  be  multiplied 
by  [(1  +  cos2  20g)  /  (1  +  cos  26b)]1'  . 

Several  ccaiments  should  be  made  at  this  point.  The  strain  and  damage 
distributions  are  expected  to  be  continuous  and  this  can  be  approximated 
by  subdividing  the  ion- implanted  layer  into  a  large  number,  N,  of  laminae. 
However,  due  to  the  relatively  large  breadth  of  the  Bragg  peaks,  (about 
0.1  degrees),  it  is  found  that  finer  division  than  about  150A  per  lamina 
does  not  appreciably  change  the  calculated  intensity  distribution.  Thus, 
a  1500A- thick,  ion- implanted  layer  can  be  well  represented  by  10  or  fewer 
discrete  laminae.  By  the  use  of  a  computer,  the  correct  strain  and  damage 
distributions  may  be  synthesized  by  trial  and  error  until  the  intensity 
matches  the  measured  rocking  curve.  Analytical  approaches,  such  as 
Fourier -analyzing  the  rocking  curve,  are  more  elegant,  but  also  more  limi¬ 
ted.  By  calculating  the  Patterson  series  of  the  rocking  curves  of  diffused 
and  ion- implanted  layers,  Afanasev,  et  al  (Ref  16)  obtained  the  average 
strain  and  thickness,  but  not  the  strain  and  damage  distributions. 

If  the  strain  and  damage  distributions  are  to  be  synthesized  by  the 
fitting  of  experimental  data,  the  sensitivity  of  the  rocking  curve  to 
small  changes  in  the  strain  and  damage  distributions  should  be  examined. 

On  theoretical  grounds  the  rocking  curve  should  be  very  sensitive  to  such 
changes.  For  a  typical  total  thickness  of  several  thousand  A,  the  relative 
phase  of  amplitudes  originating  from  differently  strained  laminae  can  vary 
by  several  times  rr.  A  slight  change  in  the  strain  distribution  should 
result  in  drastic  changes  in  the  rocking  curve.  This  is  demonstrated  in 
Figure  2(a)  where  a  rocking  curve  (solid  line)  corresponding  to  the  strain 
(solid  line)  in  Figure  2(b)  is  plotted.  The  rocking  curve  as  well  as  the 
strain  and  damage  (dotted  line)  distributions  are  typical  of  200KeV  Ne+ 
implants  in  garnet.  If  the  strain  distribution  is  changed  by  a  small 
amount,  yielding  the  dashed  distribution  in  Figure  2(b) ,  the  result  is  a 
marked  redistribution  in  the  oscillations,  as  shown  by  the  dashed  line 
in  Figure  2(a).  The  positions  and  the  relative  heights  of  the  maxima 
have  changed,  although  the  total  area  under  the  curve  has  remained 
constant.  A  similar,  but  less  pronounced,  effect  is  obtained  if  the 
damage  distribution  is  varied.  These  considerations  indicate  that  a  good 
fit  of  an  experimental  curve  yields  a  strain  distribution  which  is  accurate 
to  better  than  5  percent  at  any  depth.  Since  the  rocking  curve  is  less 
sensitive  to  damage,  the  accuracy  for  the  damage  distribution  is  only  15 
percent  or  better. 

The  synthesis  approach  becomes  very  useful  when  rocking  curves  are 
taken  after  progressively  deeper  etching,  performed  in  steps  of  100  to 
200a.  Starting  with  the  deepest  etch,  the  distribution  in  this  portion 
of  the  ion-implanted  layer  is  easily  determined  by  one  or  two  discrete 
laminae  whose  total  thickness  is  fixed.  The  strain  and  damage  in  the  next 
lamina  are  found  by  fitting  the  appropriate  rocking  curve  without  changing 
the  distribution  obtained  at  the  previous  step.  In  this  manner,  the 
distribution  is  built  up  from  the  bottom.  The  assumption  that  a  given 


lamina  is  not  appreciably  affected  by  the  presence  of  adjacent  laminae 
is  supported  by  the  agreement  between  calculated  and  measured  rocking 
curves . 

After  sane  experimentation,  hot  concentrated  nitric  acid  was  chosen 
as  a  suitable  etchant  for  implanted  garnet.  At  90#C,  implanted  material 
etch  rates  are  in  the  neighborhood  of  100  to  200A/min,  which  provides  a 
convenient  and  reproducible  means  for  removing  implanted  garnet.  The 
surface  roughening  caused  by  etching  under  these  conditions  is  within 
tolerable  limits.  The  thickness  of  implanted  material  removed  is  deter¬ 
mined  from  the  unimplanted  etch  rate  (about  25A/min  at  90°C) ,  plus  a 
mechanical  measurement  of  the  step  height  which  develops  at  the  unimplanted - 
implanted  boundary. 
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SECTION  III 


RESULTS  AND  DISCUSSION 


This  section  describes  the  key  properties  of  bubble  materials 
implanted  with  single  energy  Ne+,  He*  or  H+,  as  well  as  double  and 
triple  energy  Ne+  and  He+  combinations.  These  species  were  selected 
because  of  their  widespread  use  for  hard  bubble  suppression  and  contiguous 
disk  device  fabrication  (Refs  3,12,18).  The  device  requirement  of  low 
level  of  damage,  uniformly-strained  drive  layers  (Ref  3)  focused  the 
investigation  on  the  dependence  of  strain  and  damage  distributions  on 
implant  species,  energy,  and  dose.  In  addition,  FMR  measurements  of  some 
of  the  samples  yielded  values  for  HK,  y,  and  a. 

The  bubble  films  used  were  YSmLuCaGeIG  with  thicknesses  varying 
between  1  and  3um.  Implantation  was  performed  at  a  lpA/cnr  dosage  rate 
in  order  to  minimize  self -annealing  effects.  Multiple  energy  implanta¬ 
tions  were  made  in  order  of  descending  energy.  The  temperature  of  the 
samples  during  implantation  was  estimated  to  be  between  150  and  200°C. 

Figure  3  shows  the  (444)  CuKq  rocking  curves  of  samples  implanted 
with  200  KeV  120  KeV,  200  KeV  ♦  SO  KeV,  and  200  KeV  +  120  KeV  +  50  KeV, 
6xlO^Ne+/cnr  for  each  energy.  The  curves  corresponding  to  single-energy 
implantations  are  qualitatively  the  same,  except  that  the  200  KeV  curve 
is  sharper  and  more  intense.  The  bottom  two  curves  are  qualitatively 
different  from  each  other  as  well  as  from  the  single-energy  curves.  The 
kinematical  model  interpretation  of  these  and  other  curves  in  terms  of 
strain  and  damage  distributions  is  shown  next. 

Figure  4(a)  shows  the  experimental  and  theoretical  rocking  curves 
of  the  120  KeV,  6xlO^Ne+/cm‘'  single  implantation.  The  strain  and  damage 
distributions  which  yield  the  theoretical  rocking  curves  are  shown  in 
Figure  4(b).  Tlie  vertical  dotted  lines  in  this  figure  indicate  the  etch 
steps  at  which  rocking  curves  were  taken.  The  very  good  agreement  between 
the  experimental  and  calculated  rocking  curves  implies  that  the  strain  and 
damage  distributions  are  accurate  to  within  a  few  (<5)  percent. 

For  this  dosage,  the  damage  is  very  low  while  the  strain  is  fairly 
uniform.  The  maximum  strain  (0.47  percent)  is  at  about  1000a  from  the 
surface;  this  is  slightly  deeper  than  half  of  the  implanted  layer  thickness 
(1800A).  At  the  layer-bulk  interface  and  at  the  surface,  the  strain  values 
are  0.29  percent  and  0.34  percent,  respectively,  which  gives  the  distri¬ 
bution  a  slightly  skewed  appearance. 

The  rocking  curves  and  the  strain  and  damage  distributions  of  the 
200  KeV,  6xloJ^Ne+/cm^  sample  are  shown  in  Figures  5(a)  and  (b),  respect ivcl 
Qualitatively  the  distributions  are  the  same  as  for  the  lower  energy  case. 
The  important  differences  are  in  the  layer  thickness  (3100A)  and  maximum 
strain  (0.42  percent).  It  appears  that  the  120  KeV  implantation  results 
in  a  slightly  higher  strain  and  damage  than  the  200  KeV  implantation,  even 
though  the  dosages  are  the  same.  This  is  attributed  to  increased  straggling 
by  higher  energy  implantation  ions. 
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4  200  KeV  ♦  120  KeV  ♦  50  KeV 
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Angle  Referred  to  LPE  Peak  (deg.) 


Figure  3.  Rocking  Curves  Resulting  from  Different  Single 
and  Multiple  Inplantations. 


Theory 


Although  the  damage  is  quite  low  and  the  strain  is  uniform  in  the 
cases  discussed  above,  the  value  of  the  strain  is  barely  sufficient  to 
overcome  the  growth- induced  anisotropy  energy,  Ky.  Hard-bubble  suppression 
measurements  indicate  that  a  substantial  increase  in  the  strain  is  desirable. 
A  three- fold  increase  is  obtained  by  the  200  KeV,  ZxlO-^NeVcm-  implantation 
shown  in  Figures  6(a)  and  (b) .  The  maximum  strain  here  is  1.3  percent, 
but  its  distribution  is  much  less  uniform,  while  the  damage  has  increased 
by  much  more  than  a  factor  of  three.  Ferromagnetic  resonance  measurements 
indicate  that  for  this  damage  the  magnetic  properties  of  the  ion- implanted 
layer  are  essentially  unchanged:  4mM  =  297G, y=  (1.75+0.15)xl07/sec-0e,  and 
a  =  0.135+0.07  for  both  the  bulk  LPE  crystal  and  the  implanted  layer. 

The  bulk  H.  =  1900±50  Oe  (K  =  22.5  Kerg/cm2) ,  while  the  corresponding 
value  of  Hr  for  the  implanted  layer  is  -1090+50  Oe.  The  implantation- 
induced  anisotropy,  K,  =-35.4  Kerg/cm^.  Thus,  the  ratio  K^/K^  =  -1.6. 

The  relation  between  and  the  strain: 


3  E 
1  “  1 


1+v 


d-dn 
X  0 

Am  nr 

0 


(6) 


yields  a  value  of  -1.2x10  for  the  magnetostriction  constant,  L,,. 

Here  Young's  modulus,  E,  and  Poisson’s  ratio,  v,  were  assumed  to11 
be  2xl012  and  0.29,  respectively.  It  should  be  noted  that  while  E  and  v 
are  constant  for  most  bubble  films,  X......  is  composition-dependent  and  can 

vary  by  a  factor  of  2  or  3  in  either  airection. 


Without  an  increase  in  X-...,  the  value  of  | K- 1  cannot  be  raised 
significantly  by  additional  implantation.  A  factor  of  two  increase  in  the 
ion  dosage,  i.e.,  4xl014Ne  /cur,  becomes  self-defeating  since  the  damage 
rapidly  approaches  100  percent,  especially  in  the  middle  section  of  the 
implanted  layer.  Jhis  is  associated  with  the  creation  of  a  nonmagnetic 
region,  about  1500A  wide  (Ref  19),  which  does  not  contribute  to  the  flux- 
keeper  effect  or  to  the  charged  wall  in  contiguous  disk  devices.  Further 
increases  in  the  ion  dose  result  in  an  entirely  nonmagnetic  layer.  The 
maximum  strain  with  tolerable  damage  for  Ne+  implantations  is  about  1.5 
percent . 


We  return  now  to  the  profiles  obtained  by  multiple-energy  implantations . 
The  strain  and  damage  distributions  resulting  from  the  200  KeV  +  50  KeV, 
6xl01^Ne+/cm2  for  each  energy,  implantations  are  shown  in  Figure  7(b), 
together  with  the  appropriate  rocking  curves  in  Figure  7(a).  In  the  first 
2000A  from  the  interface  with  the  bulk,  the  strain  and  damage  are  identical 
to  those  obtained  by  the  single  200  KeV,  6xl01-5Ne+/cnr  implantation.  This 
portion  of  the  implanted  layer  is  unaffected  by  the  50  KeV0 implantation. 

The  strain  and  damage  are  increased  only  in  the  first  1100A  from  the  surface, 
and  the  shape  of  the  distribution  in  this  region  is  qualitatively  the  same 
as  for  a  single  50  KeV  implantation  (except  for  the  presence  of  a  constant 
positive  offset) .  For  this  dosage,  the  contributions  of  the  50  KeV  and 
200  KeV  implantations  are  additive,  yielding  a  maximum  strain  of  0.67 
percent  and  a  maximum  damage  of  11  percent. 
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The  rocking  curves  of  the  triple -implanted  sample  are  shown  in 
Figure  8(a).  Since  the  strain  and  damage  distributions  for  this  case 
were  expected  to  be  quite  different  from  the  single  and  double  implantations, 
rocking  curves  were  taken  at  shorter  etch  intervals  to  facilitate  the  fitting 
procedure.  Here  the  agreement  between  experiment  and  theory  is  not  as 
good  as  for  the  previous  cases ,  consequently  there  is  some  uncertainty 
(+10  percent)  in  the  distribution  in  the  first  1000A  from  the  surface. 

An  additional  rocking  curve  taken  at  an  etch  step  of  50 0.A  from  the  surface 
would  have  aided  in  obtaining  a  better  fit. 


The  strain  and  damage  distributions  shown  in  Figure  8(b)  increase 
almost  linearly  with  distance  from  the  interface  up  to  about  2500.A.  The 
maximum  strain  and  damage,  1  percent  and  31  percent,  respectively, 
indicate  that  the  contribution  of  each  implantation  is  additive. 


The  properties  of  Ne+-implanted  layers  can  be  summarized  at  this  point. 
The  layer  thickness  increases  roughly  linearly  with  ion  energv,  ranging 
from  lOOOA  t<^  3200A  for  50  KeV  and  200  KeV,  respectively.  Up  to  doses  of 
4xlO^Ne+/cm“ ,  the  maximum  strain  is  linear  with  dosage,  hut  the  strain 
distribution  becomes  progressively  nonuniform.  Relatively  uniform  distri¬ 
butions  are  qbtained  by  double-energy  implantations,  such  as  200  KeV, 
2xl0^N'e+/an“  +  50  KeV,  6x10^ JNe+/cm  ,  which  yield  a  strain  of  1.3  percent 
throughout  most  of  the  layer.  The  limiting  factor  on  the  useful  maximum 
strain  (1.5  percent)  is  the  damage  which  increases  faster  than  linearly 
with  dosage.  For  strains  less  than  1.5  percent,  the  -JrM  ,  y ,  a,  and  the 
coercivitv,  H . ,  are  relatively  unaffected  bv  implantation,  a 1  though  the 
exchange  constant,  A,  is  reduced  (Ref  14).  Depending  on  the  magneto¬ 
striction  constant,  \j^>  ratios  of  1.5  for  K- /K  are  obtainable. 


Contiguous  disk  devices  require  an  implanted  (driver)  Jayii  t  hi. -’n r  <s  , 


to  an  unimplantcd  (storage)  layer  thickness,  h 


ra  1 i  o 


tc  0.5 


(Ref  3).  For  use  with  a  storage- layer  thickness  about  !.5;.m,  a  thicker 
driver  layer  is  required  than  200  KeV  Ne+  can  produce.  For  dcepei  pene¬ 
tration,  lighter  species,  such  as  He+,  H+,  or  lit  are  more  appropriate, 
although  the  dosage  required  for  the  same  strain  can  be  two  orders  of 
magnitude  higher.  The  strain  and  damage  distributions  for  these  species 
are  not  the  same  as  for  Ne+,  and  their  effect  on  magnetic  properties  may 
be  different. 


The  rocking  curves,  and  the  strain  and  damage  distributions  for  a  60 
KeV,  5.5xl0l%levcrrr  implantation  arc  shown  in  Figures  9(a)  and  (b) .  The 
fit  for  curve  0  is  quite  good,  but  there  is  a  slight  phase  discrepancy 
for  the  remaining  curves.  This  discrepancy  corresponds  to  a  5  percent 
decrease  in  the  strain  distribution  after  the  first  etch  step.  Die  cause 
of  this  is  not  clear;  possibly  some  annealing  occurred  during  etching  at 
90°C  (for  10  min.)  or  the  removal  of  the  top  700A  relieved  sane  of  the 
strain.  This  effect  is  occasionally  observed  in  other  samples,  including 
Ne+- implanted  ones,  hut  its  occurrence  is  erratic.  In  any  case,  the  strain 
distribution  shown  is  accurate  to  5  percent  at  any  depth. 
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— •  Experiment 
—  Theory 


For  this  dosage  and  energy  the  layer  thickness  and  maximum  strain  are 
comparable  to  the  200  KeV,  2xl(r4Ne+  implantation.  The  strain  distribution 
is  slightly  more  skewed  and  less  uniform.  Although  there  is  a  ±10  percent 
uncertainty  in  the  damage  distribution,  its  peak  value  is  about  50  percent 
less  than  for  Ne+.  This  result  is  important  to  the  physics  of  ion-lattice 
interaction,  but  within  the  resolution  of  our  measurements,  does  not  affect 
the  garnet's  magnetic  properties  differently.  The  values  for  4ttM  ,  y,  a, 
and  are  essentially  the  same  as  obtained  for  Ne+.  " 

The  penetration  depth  of  200  KeV  He+  ions  in  garnet  is  0.75um,  as 
shown  in  Figure  10(b).  The  dosage  is  lxlO-^HeVcm^ ,  resulting  in  1.6 
percent  maximum  strain  and  46  percent  maximum  damage.  Note  that  the 
distance  scale  has  been  increased  by  a  factor  of  2,  so  that  the  distribution 
is  not  significantly  more  peaked  than  that  shown  in  Figure  9(b).  The  rocking 
curve  in  Figure  10(a)  is  theoretical,  but  the  agreement  with  experiment  is 
reasonably  good.  The  uncertainty  in  the  distributions  is  not  more  than  +10 
percent.  It  is  clear  from  the  strain  distribution  that  multiple  energy 
implantations  are  necessary  in  order  to  obtain  a  uniformly  strained  layer. 

Ion -implantation  depth  is  roughly  proportional  to  the  ion-accelerating 
voltage  and  to  the  reciprocal  of  the  mas|  of  the  implanted  ion.  Thus,  the 
implanted  layer  thickness  for  200  KeV  Ne  ,  He+,  and  H+  was  found  to  be  0.32, 
0.75,  and  1.25ym,  respectively.  The  implantation- induced  strain  and  damage 
are  approximately  symmetrically  distributed  through  the  implanted  layer 
for  Ne+  but  are  skewed  toward  the  end  of  the  penetration  range  for  lighter 
ions. 

Despite  the  need  for  deeper  implantations,  studies  on  H+-implanted 
garnets  have  not  been  very  numerous  because  of  the  high  dosage  needed,  and 
the  lack  of  a  suitable  etchant.  For  a  strain  of  about  0.4  percent,  the 
dosage  requirements  at  200  KeV  are  6xlO^Ne+/cm  ,  1.5x10^1  leY/cm2  and 
2xl(U,H+/cirr .  The  durations  of  implantation  for  He+  and  H+  are  therefore 
about  25-  and  300 -times  that  for  Ne  at  the  same  ion-beam  current.  Higher 
ion-beam  currents  may  be  employed  to  reduce  the  light-ion  implantation  time 
(within  equipment  limitations) ,  but  an  accompanying  rise  in  wafer  surface 
temperature  produces  self -annealing  which  is  counterproductive  to  increasing 
strain.  In  order  to  analyze  the  strain  and  damage  distributions  by  the 
techniques  developed  for  Ne+- implantations,  a  uniformly -acting  slow  etchant 
is  required.  Unfortunately,  hot  HN0-  attacks  H+-implanted  material  unevenly, 
producing  a  roughened  surface  so  that  successive  rocking  curves  cannot  be 
prepared.  It  is  possible  to  compare  the  as -implanted  rocking  curves  of  Ne+- 
and  H+-implanted  garnets  to  obtain  the  maximum  strain,  and  estimates  of  the 
strain  and  damage  distributions  in  the  H^-implanted  material.  Detailed 
measurements  however  must  await  the  development  of  a  suitable  etching 
technique . 

Our  results  indicate  that  uniformly  strained  layers  with  appropriate 
thickness,  magnetic  properties  and  low  damage  can  be  produced  by  double-  or 
triple-energy  implantations  of  Ne+  and  He+.  Since  the  dosage  requirements 
for  H+  exceed  those  for  He+  by  five-fold,  H+  implantation  is  less  attractive. 


SECTION  IV 


CONCLUSIONS 


A  method  for  using  X-ray  diffraction  to  study  ion- implanted  crystals 
has  been  described.  Strain  and  damage  distributions  in  bubble  garnets 
for  single  energy,  and  in  some  cases  multiple  energy,  Ne+,  He+,  and  H+ 
implantations  were  discussed.  Sufficiently  uniform  strains  of  appropriate 
values  (1.0  to  1.5  percent)  are  obtainable  by  multiple  energy,  dosage, 
and/or  species  implantation.  The  maximum  usable  strain  is  limited  by 
the  damage,  which  increases  rapidly  for  strains  over  about  1.5  percent. 

For  low  (<35  percent)  damage,  the  magnetic  properties  of  the  garnet  are 
essentially  unchanged,  except  for  the  creation  of  the  stress -induced 
anisotropy  energy,  K. .  This  parameter  is  proportional  to  the  product  of 
the  strain  and  the  magnetostriction  constant,  L...  For  contiguous  disk 
device  applications,  the  value  of  K.  should  exceea  the  growth -induced 
anisotropy  K  by  ~50  percent.  Ferromagnetic  resonance  measurements  on 
1. 5pm -diameter  bubble  materials  indicate  values  of  about  50  Kerg/arrr  and 
about  75  Kerg/cnm  for  the  and  K^,  respectively. 

As  the  bubble  size  is  decreased  for  higher  bit  densities,  both  4nM 
and  K  increase,  making  it  more  difficult  to  create  a  K.  of  sufficient  s 
magnitude.  The  composition-dependent  X...  could  be  raided  to  meet  this 
problem,  but  compositional  changes  affect  other  magnetic  parameters 
including  K  .  Thus  for  smaller  bubble  material,  it  may  be  necessary  to 
resort  to  double -layered  films  that  have  a  driver  layer  with  a  low  K 
(which  can  be  exceeded  by  K.) .  However,  fabrication  costs  for  double- 
layer  films  will  be  higher  than  for  single-layer  films,  and  the  possibility 
of  defect  formation  will  increase.  Another  aspect  of  smaller-bubble  material 
is  that  the  required  thickness  of  the  implanted  layer  also  decreases,  making 
a  single-  or  double-Ne+  implantation  suitable. 

A  more  complete  understanding  of  implanted  layer  properties  should 
result  from  the  use  of  additional  observation  techniques.  Ferrofluid 
decoration  would  permit  estimation  of  the  strength  and  length  of  the 
charged  wall,  and  this  information  could  be  correlated  to  the  driver  layer 
thickness,  strain,  and  damage.  In  addition,  independent  determination  of 
Xm  by  FMR  on  unimplanted,  externally- stressed  samples  would  aid  in  deter¬ 
mining  the  potential  of  a  material  for  a  high  K.  prior  to  implantation. 
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